A study of random laser modes in disordered photonic crystals 
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We studied lasing modes in a disordered photonic crystal. The scaling of the lasing threshold 
with the system size depends on the strength of disorder. For sufficiently large size, the minimum 
of the lasing threshold occurs at some finite value of disorder strength. The highest random cavity 
quality factor was comparable to that of an intentionally introduced single defect. At the minimum, 
the lasing threshold showed a super-exponential decrease with the size of the system. We explain it 
through a migration of the lasing mode frequencies toward the photonic bandgap center, where the 
localization length takes the minimum value. Random lasers with exponentially low thresholds are 
predicted. 



A finite open system of scattering particles can be 
characterized by a set of quasi-stationary (leaky) optical 
modes. When gain, sufficient to compensate the leak- 
age in at least one mode, is introduced in such system - 
random laser is formed | |, |, |, |, @, |, §. Uncon- 
ventional feedback mechanism in a random laser leads to 
properties interesting from both fundamental and prac- 
tical point of view |, 1 1, |, % 0, 1 1, 0, 0, 0. 
While previous studies concentrated mainly on disor- 
dered systems, the resent experiments by Shkunov et al 
[0 demonstrated random lasing peaks associated with 
the bandgap in partially ordered system. In this let- 
ter, we obtain the lasing threshold for a random laser 
with different degree of ordering. Important information 
about these states can be extracted from the scaling of 
the average lasing threshold with the system size [ p^ , 
L. For increasing values of the disorder parameter we 
find 5 scaling regimes: (a) photonic band-edge, l/L^, (b) 
transitional super-exponential, (c) bandgap-related ex- 
ponential, (d) diffusive, and (e) disorder- induced 
exponential, due to Anderson localized modes, regimes. 
In this letter we show that by optimizing the disorderness 
of the sample one can dramatically reduce the threshold 
of random laser, to the values comparable to those of pho- 
tonic bandgap defect lasers [|l5j. This finding can open 
up a road to practical applications of random lasers. 

Finding the lasing threshold of a random laser theo- 
retically is a difficult problem. In disordered photonic 
crystals one cannot take simplifying assumptions such 
as independent scattering approximation (low density 
limit), neglect the finite size of the scatterers, or the con- 
fined dimensions of the system. Random laser modes 
have the highest Q's and therefore weakly contribute to 
the transport properties of the system. Among other 
methods finite difference time domain (FDTD) 

method has been shown to be a convenient tool in study- 
ing random laser modes in ID |0 and 2D It was 
recently shown [l^ that above the threshold, the 
shape of the lasing modes remain the same as in pas- 
sive medium. In this work, we used FDTD method to 
find the modes with the smallest energy decay rate in an 
open passive 2D random medium system with different 



degree of order. These modes determine the value of the 
lasing threshold in an active medium. After an initial 
pulsed excitation with the frequency coe in vicinity of the 
photonic bandgap, electric energy of the system followed 
£ (X Re [exp2ia'm(l + i/'2Qm)] dependence. Um and Qm 
were identified as the frequency and quality factor of the 
random cavity. In what follows, we will refer to 1/Qm 
as the lasing threshold of the system. Our main results 
can be summarized as follows: (i) for sufficiently wide 
bandgaps, the minimum of the lasing threshold develops 
at some finite value of disorder parameter; (ii) at this 
"optimal" disorderness, the cavity Q is determined by 
the localization length similar to that of single defect in 
the ordered sample, leading to a similar random cavity 
Q; (iii) with an increase of the system size the optimal 
disorder decreases; (iv) near the minimum, lasing thresh- 
old scales super-exponentially on sample size, which we 
ascribe to the migration of mode frequency toward the 
bandgap center and, therefore the decrease of the "typi- 
cal" localization length. 

We consider a 2D photonic crystal made oi N cc 
cylinders with diameter d — 98nm and refractive index 
no ~ 2.2. The cylinders were arranged into hexagonal 
lattice with nearest-neighbor distance of a = lAQnm. In 
the absence of disorder, the infinite system with these pa- 
rameters has full bandgap in [361nrn, 426nm] range for 
TM modes (electric field along the cylinder axis). The 
disorder in the system was introduced in 2 ways: by 
uniformly randomizing the refractive index of different 
cylinders in the range [n^ — w„(no — 1), no -I- WnijiQ — 1)] 
and diameter - [d{\ — Wd),d{l + Wd)]- Care was taken 
in order to avoid the uncontrollable disorder due to dis- 
cretization of the grid. Disorder in the system was char- 
acterized with parameter Se =< J (e(r) — £o(r))^dr >i/2 
/ (/ £Q(r)dr)^/^, where £o(r) and e{r) is the dielectric con- 
stant distribution in ordered and disordered samples re- 
spectively. < ... > stands for the average over different 
disorder configurations. In the present work we studied 
the systems with 11 different disorder strengths: 1 to 10 
had Wn G [0.1,1.0] with the increment 0.1 and Wd = 0, 
the 11th had Wn — 1.0 and Wd — 0.43. This led to 
variations of dielectric constant, Se, from weak - 0.08 to 
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FIG. 1: Average inverse random cavity quality factor normal- 
ized by total number of scatterers A'^ oc L^, as a function of 
disorderness. 

strong - 0.95 disorder. Later we will discuss the effects of 
this particular choice of these types of the disorder. We 
used standard FDTD method to study the evolution of 
TM electromagnetic field in the described system system 
with A/20 resolution grid. In the studied range of system 
sizes, N £ [75,600], the required grid sizes varied from 
200 X 200 to 500 x 500. To mimic an open system, a 
buffer layer of air of 150nm thickness was kept around 
the sample, followed by uniaxial perfectly matched layer 
I p^ . To excite the system we initially launched a short 
~ 10/s pulse at every grid point. The frequency of the 
pulse was at the center of the bandgap, 391nm, of the or- 
dered structure. In the frequency domain the full width 
at half maximum of the seed pulse was of the order of the 
bandgap width. Thus the seed pulse excited all states 
within the stop band and near the band edges. After 
initial excitation pulse, the system was left on its own. 
Field distribution in the system and total electric energy 
£ = 1/2 / £(r)E^(r)c?r was monitored. 

Right after the initial pulse, the competition between 
the modes |l^ with the different lifetimes led to 
the complicated evolution of £. However, after a suf- 
ficient time only the mode with the longest lifetime 
(maximum quality factor) survived. £ followed a mono- 
exponential dependence Re[exp2iujmi^ + i/'^Qm)], from 
which we extracted the maximum quality factor Qm and 
frequency Um for this particular realization of the dis- 
order. At the same time the spatial pattern, E(r), sta- 
bilized and the mode profile could be seen. Generally, 
the time needed to reach the mono-exponential regime 
varied from about 0.5ps for the smallest system to lOps 
for the largest. Finally, the 1/Qm, representing lasing 
threshold, was averaged over 1000 {N — 75), or 100 
(N = 137, 188, 261, 368, 449, 608) disorder realizations. 

Fig. IT] shows the dependence of the average lasing 



FIG. 2: Scaling of the average inverse quality factor with the 
lateral size of the system L oc N^^'^ in 5 regimes: band- 
edge type - circles; super-exponential - stars (line is used to 
connect points); bandgap exponential - squares; diffusion 
type - triangles; and disorder exponential - inverted triangles. 
These correspond to «;„ =0.1 {Se ~ 0.08), Wn = 0.2 {Se ~ 
0.15), Wn = 0.3 {Se ~ 0.23), w„ = 0.7 {Se ~ 0.53), and w„ = 
1.0 &i Wd = 0.43 {Se ~ 0.95) values of disorder parameter 
respectively. 



threshold normalized by iV as a function of the disor- 
derness, different curves correspond to different system 
sizes. This particular normalization makes it easy to 
see the deviation from diffusion predicted ||] dependence 
1/L^ oc N. One can see that significantly different scal- 
ing at different Se leads to the appearance of a minimum 
at the finite disorder strength. Fig. || plots L dependence 
of the threshold for the fixed disorder strengths. 

The understanding of this behavior comes from the po- 
sition of the mode frequencies - Fig. ^. For small 5e the 
frequencies are concentrated at lower (long-wavelength) 
band-edge, and they (as well as Qm) are independent 
on the position uje of the excitation pulse. The rea- 
son for that is the way the disorder was introduced into 
the system. The long-wavelength modes are mostly con- 
centrated in the high dielectric medium, which is be- 
ing disordered by the index fluctuations. At w„ = 0.1 
the mode frequencies fell in the immediate vicinity of 
the band-edge. Fig. ||a. Lasing from the band-edge 
modes is well studied in the case of the ordered struc- 
tures [[2[ |0|, H, H, with l/L^ dependence of 
the threshold ||l2] giving a good fit with our results. At 
the increased disorder, w„ = 0.2, we observe the decrease 
of the threshold that is noticeably faster than exponen- 
tial, expected for locaHzed modes 0. This unexpected 
behavior becomes clear from Fig. Bp. The value of the 
threshold can be estimated as 1/Q oc exp(— L/2^(tt'm)), 
where ^{ojm) is some "typical" value of the localization 
length at the frequency uJm- From Fig. ^ one can see 
that with increased size, the frequency of the modes ad- 
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vance toward the band center, where ^ is the shortest, 
magnifying the exponential dependence of the threshold. 
The explanation for this frequency migration with the 
increase of the size comes from the fact that in the finite 
system the modes deep into the bandgap are not acces- 
sible due to low density of states there. The Urbach- 



2^. The support for 
J, where 1/Q follows 



like behavior can be expected 111:, 
this, indirectly comes from Fig. 
rather well the exponential dependence (solid line) on 
the frequency shift from the band-edge of the ordered 
system. However, because of the special nature of the 
modes no direct conclusions can be made. Assuming 
Urbach-like dependence of the density of states, the ad- 
vancement of the mode frequency can be estimated from 
the condition that total number of states, N, times the 
probability of having a state located Aw(iV) away from 
the band-edge, exp(— a(fc) ■ Auj{N)), is still of order of 
one. Here, a{Se) is exponential slope in the density of 
states, that should decrease with the increase of the dis- 
order Se. For small disorders a~^{6e) ^ AEpbg, where 
AEpBG is the width of the photonic gap. Therefore, 
for weak disorder (or small system size) the bandedge- 
type modes have the highest Q. The crossover to the 
super-exponential dependence of the threshold occurs 
when the cavity Q associated with the localized states 
with the smallest localization length ^{Auj{N)) avail- 
able for this size N exceeds that of the bandedge-type 
mode: eiip[-N^/^a/2^{Auj{N))] ~ Stronger 
size dependence in the super-exponential regime means 
that even A^~3/2 band-edge type dependence observed 
at smaller disorder would eventually switch to the 
super-exponential one. However, the latter can be ex- 
pected to saturate at larger sizes or disorders when 
the mode frequency reaches the bandgap center: N ■ 
exp{—a{Se) ■ AEpbg 1"^) ~ where the localization 
length is the smallest. Therefore, we expect the lim- 
iting scaling of the lasing threshold to be exponential 
e-xp[-N^/'^a/2i{AEpBG/2)]. Three important conclu- 
sions (i-iii) listed in the introduction immediately follow, 
(iii) is justified because (^{AEpbg/'^) is the smallest in 
the ordered sample. The saturated exponential depen- 
dence can be seen in our example for Wn = 0.3 - squares 
in Fig. |. 

In our system the sharp rise of the threshold at Se ~ 0.4 
on Fig. 1^ is attributed to the removal of the bandgap. 
This could be seen as the loss of the hexagonal symmetry 
of the observed modes as well a sensitivity of the modes 
to the excitation pulse position uje- In this regime, the 
frequency of the modes is not associated with the pho- 
tonic gap, which does not exist anymore. However, in 
order to make a direct comparison with the ordered case, 
we kept the excitation pulse the same as before. The 
exact cause of the disappearance of the bandgap is be- 
ing debated in the literature ^ |2^, and is not the 
subject of this study. In our particular case we found a 
simple explanation of the threshold behavior in the way 
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FIG. 3: 1/Q vs cavity mode frequency for 20 realizations 
of disorder. Circles, squares, and triangles correspond to TV 
equal to 75, 188, and 608 respectively. Four graphs correspond 
to different disorder parameters: Wn =0.1 {Se ~ 0.08) - (a), 
w„ = 0.2 (Se ~ 0.15) - (b), Wn = 0.6 {Se ~ 0.45) - (c), 
w„ = 1.0 & Wd = 0.43 {Se ~ 0.95) - (d). 



the disorder was introduced. Indeed, the fluctuating in- 
dex of refraction leads to the fluctuation of the frequency 
of the Mie resonances of the particle. For box distribu- 
tion, there exists a value of w„ = 0.6 when the defect Mie 
resonance falls into the gap. Fig. ^. This value matches 
the value of disorder parameter, when the sharp increase 
of the lasing threshold is observed on Fig. |l|. Moreover, 
Fig. ||c shows that at this crossover disorder, the modes 
avoid the region of high scattering, which is the conse- 
quence of the sharp boundary in the distribution of n. 
It also indicates the presence of the residual bandgap, 
where the lasing frequencies are concentrated. 

At Wn > 0.7 the photonic bandgap ceased to exist. 
Triangles on Fig. || correspond to i/7„ = 0.7, and were 
successfully fitted with diflFusion [Q, |^, ||, ||, |l^ scaling 
dependence Deviations from this dependence can 

be seen from Figs. |^,^ at the largest sizes studied, where 
L > £,Anderson the statcs bccomc localized again by virtue 
of Anderson localization ||, 12 . The exponential depen- 
dence becomes especially pronounced at the largest disor- 
der studied (see inverted triangles on Fig. ||), where the 
transition from to exponential dependence comes 

at small system sizes. We want to point out that even at 
such strong disorder, the obtained modes had a collective 
nature, rather than the single particle's (high-order res- 
onances), which were concentrated at higher frequencies. 
Comparing the localization length of these states to that 
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FIG. 4: Mie scattering efficiency vs wavelength for cylinders 
of diameter 98nm. Refractive index n — 2.2 - solid line. For 
Wn = 0.6 {5e ~ 0.45), the refractive index fluctuates from 
n = 1.48 (dash-dotted line) to n = 2.92 (dashed line) 



of bandgap nature we see the factor of two difference, 
which makes the latter preferable - Fig. ^ This is in 
accordance with experiments of Ref. (llj . 

In conclusion, we studied the modes responsible for 
the lasing threshold of a random laser based on disor- 
dered photonic crystal with a complete bandgap. At dif- 
ferent disorder strengths, we identified 5 scaling regimes 
of the lasing threshold. Difference in scaling rates led to 
the development of the minimum at "optimal" disorder, 
where the threshold takes minimum value. We, for the 
first time, showed the regime of super-exponential scal- 
ing, which was attributed to the migration of the lasing 
frequencies to the center of the bandgap. We also predict 
that the value of the "optimal" disorder should decrease 
with increase of the volume of the system. These results 
can be generalized to the case of the 3D disordered pho- 
tonic crystals. 
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